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Abstract. In this paper we study the Hamiltonian structure of the second Painleve hierarchy, an infi- 
nite sequence of nonlinear ordinary differential equations containing PII as its simplest equation. The n-th 
element of the hierarchy is a non linear ODE of order 2n in the independent variable z depending on n pa- 
rameters denoted by t\ , . . . , t n _i and a n . We introduce new canonical coordinates and obtain Hamiltonians 
for the z and t\, . . . , t n — 1 evolutions. We give explicit formulae for these Hamiltonians showing that they 
are polynomials in our canonical coordinates. 
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1. Introduction 

In this paper we study the Hamiltonian structure of the second Painleve hierarchy, an infinite sequence 
of nonlinear ordinary differential equations containing 

Pn : w zz = 2w 3 + z w + ct\ , 

as its simplest equation. The n-th element of the hierarchy is a non linear ODE of order 2n, depending on 
n parameters denoted by t±, . . . , and a n : 



>(n) 



- + 2w)C n [t 



i=i 



where C n is the operator defined by the recursion relation 
d 



(1) 



d 3 d 



w = zw + a n 



Co\w z 



n > 1, 



The second Painleve equation and its hierarchy appear in several applications including Helc-Shaw geom- 
etry JJj, nonlinear optics ^H] and random matrix theory [SHI El to name only a few. 

The Hamiltonian structure of the classical six Painleve equations was discovered long ago by Okamoto 
|36| . Jimbo and Miwa |2H]. In the case of n = 1, i.e. Pn, the Hamiltonian is 



4 4 2 



where 



Q = 4w, 



p =b" r 
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Using such formulation Okamoto was able to describe his initial conditions space, to characterize the action of 
the Backlund transformations found by Gambier JHJ an d Lukashevich jSJ in terms of affine Weyl groups and 
to produce immediately the so called Riccati-type classical solutions of the second Painelve equation ^3 1211 
[S]. Also several properties of the Yablonskii-Vorob'ev polynomials |41l 140] describing the rational solutions 
were proved using the Hamiltonian formulation Umemura and Watanabe used the Hamiltonian 

structure in to prove the irreducibility of Pn. 

In this paper we introduce canonical coordinates Pi, . . . , P n , Q\, . . . , Q n and a Hamiltonian function 
such that Pj™' 1 is equivalent to 

dQj _ op m&) ._ 

( ' dz ~ dp ' dz ~ dQi ' % ~ 

In particular we show that TiS n ^ is a polynomial in Pi, ... , P n , Qi, . . . , Q n and that the Hamiltonian equations 
satisfy the Painleve property. 

Observe that starting from the second element of the hierarchy, the parameters t\, . . . ,i„_i appear in 

(n) 

Pjj . The solutions w will depend on the times t\, . . . , i n _i according to the equation 
(3) (2k + l)-^-+d z (d z +2w)C k [w z ~w 2 ] =0, fc = l,...,n-l, 

which is actually the fc-th element of the mKdV hierarchy. In fact, the second Painleve hierarchy was 
discovered as self-similarity reduction of the mKdV hierarchy [3|SJE1 (details on this derivation are recalled 
in section [21 below) . 

We refer to the evolution in the times ii,...,t„_i, as time-flows. We prove that the time-flows are 
Hamiltonian and we compute the Hamiltonians Ti^\ . . . , H^i such that the system 

dQ t _ dHj n) dPj _ dH^ . = ] 
dt k dPi ' dt k dQi 

is equivalent to ®. These Hamiltonians are also polynomials in Pi, . . . , P n , Qi, . . . , Q n and again the 
Painleve property is satisfied. Amazingly, we obtain explicit formulae for HS n ^ and Hi , k = 1,. .. ,n — 1, 
in terms of Pi, . . . , P n , Q\, . . . , Q n without need of recursion. 

Our construction is based on the fact that the second Painleve hierarchy can be interpreted as monodromy 
preserving deformation equation 1 of an auxiliary linear system 

^f = ^ (n) (A; z, *!,...,*„_!)* 

where is a matrix function of A holomorphic in C* , having a simple pole at and a pole of order 2n + 2 
at infinity. The isomonodromic condition is expressed by the zero-curvature conditions 

3Ai n ) BE , s //it") 

(4a) A^-ik = \*^\ B =- 



dz d\ L ' J ' V 4 ™ 

(4b) d x M {k) - (2k + l)d th A (n) = -[M( fe U W ], 

where given any Laurent series L of A, (L) + denotes its non-negative part and the relation between M k and 
is explained in Section [3] below. 
We interpret equation l|4a|) as flow on the dual space of the following twisted loop algebra: 

0- = |x(A) = ^A l A l |A(A)a 1=( r 1 A(-A)|/ 02 , i+1 , *i = ( J j) = 
92n+1 = \X(X)= J2 X i X i \X(X)a 1 = a 1 X(-X)\ 



^The isomonodromic deformation problem for the second Painleve hierarchy with ii = • • • = t n — 1 = was derived in 1101 
and in 28 following the approach proposed in starting form the isomonodromy deformation problem given in 1141 for the 
second Painleve equation. Here we generalize the construction of 1101 to the case of generic values ti, ... , tn—l- 
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Its Lie bracket and a precise description of the corresponding loop group can be found in Section 0J The 
dual space g*_ of g_ can be identified with 

g* = |s(A)=f]s J A 4 |S l e S [(2,C), S(A)ff 1 =-(7 1 H(-A)|/fl5ft +1 , 
02»+i = I 3(A) = £ S 4 A l |Hi G «I(2, C), H(A)ai = -<r x S(-A) 1 /g* 2n+1 , 

{ 2n+l ) 

by the killing form in the loop algebra sl(2, C) (see l|23|) below). 

More precisely, to interpret l|4a|) as flow on the gl, we denote by A the dynamical part of A^- n \ i.e. 
A = (U (rl) ) + , and we define 



B 



A^X 1 



4"- 

where given any Laurent series L of A, (£)- denotes its strictly negative part. In this way A G gL and 
B G g_. Then we show that equation l|4a|l is equivalent to 

,„, dA dA rry 

(5) — -— = [B,A]=bA b A, 

az az 

so that the r.h.s. defines a vector field on the coadjoint orbit Oa of A obtained by fixing the values of the 
parameters t\, . . . , t n -\ which are the Casimirs of the standard Poisson bracket on We then prove that 
the vector field defined by &d* B A, is Hamiltonian with Hamiltonian given by 

#W := - Jj-TrRes ^A 1_2n (U (n) ) 2 

where Res denotes the formal residue at 0, i.e. the coefficient of the term in A -1 . 

Let us now describe our canonical coordinates Pi, ... , P n , Q%, . . . , Q n on the coadjoint orbit Oa- Our 
construction is based on the so-called algebro geometric Darboux coordinates. The latter are the projec- 
tions of the points in the divisor of a line-bundle on the spectral curve |131 1351 0] . However this classical 
construction in our case is not straightforward: on one side there are too many points in the divisor, on 
the other side the dependence of the matrix entries of A^ on w and its z-derivatives is very complicated 
as it involves the recursive relation In particular, this means that the matrix entries are dependent on 
each other in a complicated way. We have resolved both these problems by expressing the Lenard recursion 
operator C n in terms of Cq, . . . , £ n _i and their derivatives (see (|54|l below) to obtain the matrix entries as 
polynomials of the canonical coordinates (see Theorem 16. II below). It is worth noting that in ^5] and |16|. 
a different set of canonical coordinates on the coadjoint orbits of s/(2,C) was found. Although these pre- 
serve the Painlcvc property it is difficult to apply their construction to our case because of the complicated 
dependencies between the matrix entries of A^ n ' . 

Resuming we prove that the map 

(6) A ^(Pl,...,P n ,Ql,...,Qn) 

gives a system of rational Darboux coordinates on the coadjoint orbit. However this map depends on 
z,ti,..., t n -i explicitly. This produces a shift 5H^ in the Hamiltonian H^ n \ so that the Hamiltonian HS n ^> 
appearing in J2| is given by 

n (n) = H (n) + SH < ~ n K 

Thanks to the fact that all formulae are explicit, we can compute this shift (see H66|l below). 

The idea of interpreting the isomonodromic deformation equations as Hamiltonian flows on the coadjoint 
orbits of the loop group SX(2,C) was already used by Harnad and Routhier to study the Hamiltonian 
structure of the six classical Painleve equations (see also |22])- Later Krichever used the Lax represen- 
tation approach to construct the isomonodromy equations for meromorphic connections with irregular and 
regular singularities on algebraic curves. 

On the other hand, the coadjoint orbits of the dual loop algebras can also be thought of as moduli spaces of 
meromorphic connections on Riemann surfaces. Audin generalized the Aityah-Bott symplectic structure 
on the moduli space of holomorphic connections to the case where the Riemann-surfaces can have boundaries. 
This symplectic structure was used by Hitchin to study the Schlesinger system j2Sj- Later, Boalch further 
generalized this symplectic structure to the moduli space of generic meromorphic connections. He then 
showed that these moduli spaces are isomorphic to the coadjoint orbits as symplectic manifolds, and that 
isomonodromic flows induce symplectomorphism between coadjoint orbits at different times. This abstractly 
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indicates that one could express a general isomonodromic deformation as a time-dependent Hamiltonian flow 
on the dual loop algebra. Woodhouse 021 showed that the isomonodromic deformations with poles fixed are 
autonomous Hamiltonian systems w.r.t. the Konstant Kirillov Poisson structure on a central extension of 
the dual loop algebra (which is different from what we are considering in this paper). 

Our construction is very explicit and allows us to go further to interpret the time-flows l|4bj) as flows on 
the dual loop algebra gL. This is not trivial because the term d\M^ is not tangent to the coadjoint orbit 
O a- To overcome this difficulty we introduce new coordinates u\, . . . , ui n on the coajoint orbit O a and new 
times Si, . . . , s„_i such that 

d Sk A (n) = [L k ,A {n) ]+d x L k , L k = (A {n) \ 1 - 2k ) + , d u Sk A [n) = d x L k , k = l,...,n. 

where <9" fc denotes the partial derivative with respect to s k when u is fixed. 
We then interpret the equation 

(7) d 8k AW - d x L k = (d Sk - dl)A^ = [L k ,A n \ 

(n) (n) 

as flow on the coadjoint orbit Oa arid compute the corresponding Hamiltonians h\ , . . . , h n _ l 

4") = I T rRes (V- 2fc (A^y) , k = l,...,n, 

so that in particular — — %r- Finally we show that the time-flows Hamiltonians 7"^ n \ . . . ,7i^\ are 

given in terms of hi, ■ ■ ■ , /i^i (and their shifts due to the explicit dependence of JBJ on z, ti , . . . , by 
a simple formula (see Corollary 18 . 51 below) . 

This result gives an insight into how one could express a general isomonodromic deformation as a non- 
autonomous Hamiltonian system. 

Remark 1.1. Note that h\ , . . . , h[^ are spectral invariants. In the context of iso-spectral deformations, 
this can be used to show that the algebro-geometric Darboux coordinates are separated for the isospectral 
system (see, for example |37| ). However, in the isomonodromic case, all the functions hi , . . . , hffl are 
non-autonomous, i.e. they involve the variables ti,...,f n _i and z explicitly. Therefore we don't have 
separability in the sense of classical mechanics. 

This paper is organized as follows. In Scction[31 we recall the derivation of the second Painleve hierarchy 
as self-similarity reduction of the mKdV hierarchy. In Section we describe the monodromy problem 
associated to the second Painleve hierarchy. In Section we introduce our twisted loop algebra, interpret 
equation l|4a|l as flow on its dual space gl study the Poisson bracket on 2-- We prove that the parameters 
ti, . . . , t n ^i belong to the kernel of such Poisson bracket and characterize the symplectic leaves. In Section 
|S| we introduce our coordinates Pi, . . . , P n , Qi, . . . , Q n and prove that they are canonical with respect to the 
symplectic form on the coadjoint orbit. In Sectional we obtain formula H54J1 expressing the Lenard recursion 
operator C n in terms of Co, . . . , £ n -i and their derivatives and we give the explicit formulae for the matrix 
entries of A^ n ' in terms of our canonical coordinates. We give an explicit example to illustrate our procedure 
in detail. In Section[7| we compute the Hamiltonians Ti}- n \ In Section|Sl we interpret equation @ as flow on 
the coadjoint orbit, compute the corresponding Hamiltonians h\ , . . . , n^_i and show that they are spectral 

invariants. Finally we obtain the Hamiltonians TL^, . . . , We follow all details of our construction in 

an explicit example. 

Acknowledgments. The authors are grateful to H. Flaschka, M. Talon, J. Harnad, A. Hone and to V. 
Kuznetsov for helpful conversations. This research was sponsored by the ESF grant MISGAM, the Marie 
Curie network ENIGMA and by EPSRC. 

2. The second Pn hierarchy 

The second Painleve hierarchy is obtained as self- similarity reduction of the modified Korteweg-de Vries 
(mKdV) hierarchy (see 0EHH1): 

(8) J-v + | (A + ^ Un [v x -<?] = 0, n = 0, 1,2,3, .. . 

where lZ n satisfies the Lenard recursion relation |SU] 

d ( d 3 d \ 

(9) —n n+1 =( — + A(v x -v 2 )— + 2{v x ~v 2 ) x )n n , no[u} = \. 
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Each equation in the mKdV hierarchy defines a Hamiltonian flow and can be viewed as a symmetry for all 
others. Consider the space of stationary solutions w.r.t. the symmetry defined by the n-th mKdV equation, 
i.e. the space of solutions v(x;Ti,T2, ■■■) such that g ^ v - = 0. Due to the fact that all Hamiltonian flows 
commute, all other elements of the hierarchy can be restricted to this space. 

There are also other symmetries acting on the mKdV hierarchy. They are called Virasoro symmetries. 
The n-th Virasoro symmetry is given by the following infinitesimal generator 

d :=£(2I + 1)U +1: 



ds„ f-^ dTi + i 



The stationary solutions w.r.t. this generator are by definition such that = 0. They satisfy 

£--£<a + im»£(! + *)M*->']-o, 

1=0 



and after integration 



(10) - 5^(21 + l)T l+1 ( — + 2vj1Z l [v x - v 2 } = a n , 

i=o \ x / 

where a n is some constant. 2 From the n = equation of the mKdV hierarchy, we can set T\ — —x so that 
IjlUI) is an ODE in the variable x depending on some extra parameters T 2 , . . . , T n+ i. The parameter T n+ i 
can be absorbed by the following symmetry reduction (see [5] for details): 

„ > w(z) X 

(11a) v(x,T n+1 )- 



\{2n + l)T n+1 f^ ' [(2n + l)T n+1 f^ ' 



(lib) TZi [v x - v 2 ] = — — — — — Ci[w z - w 2 }, 

[(2n + l)T n+1 } 2l/(2n+1) 

\ * t ( 2 * + 1 ) r '+ 1 ; i * 1 

(11C) t0 = ~ Z > U ^[(2n + l)T n+1 f l+ ^ + ^ l = tn = 1 - 

In this way we obtain the Second Painleve Hierarchy: 3 

(12) P|" } : ^ + 2v?j C n [w z - w 2 ] + 2w ) £ i K " ™ 2 ] = «w + On, n > 1, 

where a n are constants and C n is the operator defined by 

d / d 3 d \ 

(13) —C n+1 =[^+A(w z ~w 2 )— + 2{w z -w 2 ) z \C n , Co[w z -w 2 ] = \. 

Example 2.1. For n = 1, equation (|12f> is Pn: 

w zz — 2w 3 = zw + a.\. 

For n = 2, it is: 

ti(w zz — 2iw ) + (w zzzz — 10i/w z — 10iu w zz + 6w ) = zw + a 2 . 

In the case of the Virasoro symmetries, only the first n flows of the mKdV hierarchy can be restricted to 
the space of such stationary solutions. These give the flows in ti, . . . , t n -i: 

(14) ( 2 k + l)— + d z (d z +2w)£ k [w z ~w 2 ] =0, fc = l,...,n-l. 

dt k 

We shall call the flows in t±, . . . , t n -\ time-flows. 



2 Recently S. Kakei I2(il . proposed a new way to obtain the second Painleve equation directly as a reduction of mKdV, 
without integration. In his difference— operator formulation, the constant ct\ appears as a parameter in the symmetry reduction, 
and not as integration constant. It would be interesting to see whether this construction can be used to produce the whole Pn 
hierarchy. 

3 To obtain exact solutions of the n-th mKdV equation, one fixes the values of ti,... ,i n _i. As a consequence, often in the 
literature the second Painleve hierarchy is presented with t\ = ■ ■ ■ = t n —l = 0. We will leave ti, . . . , t n -i free to vary instead. 
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Remark 2.2. Another Painleve hierarchy containing Pn as its first element has been introduced by Gordoa, 
Joshi and Pickering |2U| by generalizing the isomonodromic deformation equations by Jimbo Miwa. However 
it is not clear if their hierarchy is different or not from the one studied in this paper. We know that 
Koike from Kyoto University is building the Hamiltonian structure of Gordoa, Joshi and Pickering hierarchy 
as confluence procedure of the Gamier systems. Once he is successful, we may try to see whether our 
Hamiltonian system and Koike 's are related by a canonical transformation. If not, the question of whether 
the Pn hierarchy considered in this paper may or may not arise as confluence limit of the Gamier system 
remains open. 



3. Isomonodromic Problem for the Pn Hierarchy 

The isomonodromic deformation problem for the second Painleve hierarchy with t\ = ■ ■ ■ = t n —\ = was 
derived in [TU| and in j^Hl following the approach proposed in starting form the isomonodromy deformation 
problem given in ^3] f° r the second Painleve equation. Here we generalize the construction of ^H] to the 
case of generic values ij., . . . , t n -x (for details see the Appendix lA")l . 

The isomonodromic deformation problem for the Pn Hierarchy is the following: 



(15a) 

(15b) 

(15c) 
where 

with 
(16a) 



a* 

dz 

a* 
Ux 

dm 



= Bt> = 
= A {n) m 



—A w 
w X 



(2k + = (M«- 



— Xz —a n 
—ot n Xz 

(d z 
(0, + 2w)£ k 



1=1 



- 2w)C k 




A. 



21 + 1 



(i6b) = 



= 4* 



(16c) 

(16d) 
(16e) 



B 



(0 

2k+l 



B 



(') 



4 fc+i 

4 fe+i d 

2 dz Vcl : 
cl 



A 2k =0, Vfc = 0,...,i, 

Q-k [w z - W 2 ] - ^- (-^- + 2w ) Cl-k-l [w z 



k = 0,...,l-l, 



2 k 



a 



m 

2k+l 



(16f) 



(17) 



— + 2wj d- k [w z - w 2 ] , 

-B 2 k+i, k = 0,...,l-l, 
B 2k , k = 0,...,l. 
The compatibility between l|15a|) and (|15b|l gives 

dA^ dB 



2w ) Hi-k-\ [w z - w 2 ] , fc = 0, 1, 
k = 1, . . . , t, 



a 



in 

2 k 



dz dX 

which gives (|12(l (see JD]). Equation (|14J) is obtained as compatibility between (|15a|) and (|15c(l : 

(18) (1 + 2k)dt k B- d z M {k) = [6,M«], 
where for brevity we put 

Finally the compatibility between l|15bll and 115c|) is 

(19) d\M {k) - (2k + l)d tk A (n) = ~[M {k \A (n) ] 



A 



The proof of the fact that equations (|17fl , (|18fl and l|19fl are indeed consistent is sketched in the Appendix 
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To simplify our computations, it is convenient to introduce some new notations. We define 

z, 



4fe+i — yi^2fc+i' k—l,...,n, a[ ) — ^ UA[ ^ 



1=1 1=1 

n 

n 

42' = X! i/B 2fc » fc = 1, . . . ,n, b ( Q ] = -a n , 
i=i 

where t n = 1, so that we can write 

,00 \2fe V^™ I,00\2fc-1 , \2fe 



(20) 



j(n) ._ ( l^k=a a 2k+i A 2^k=o°2k A +l^k=a°2k+i A 

' \ X^ n J,00\2fc-1 v»n-l ,(ra) x 2fc v-»" „00 \2fc 

\2^fc=0 O 2fc A ~ 2^fc=0 °2fc+l A ~ 2^/c=0 a 2fc+l A 

4. COADJOINT ORBIT INTERPRETATION 



In this section we show that the Hamiltonian structure of the second Painleve hierarchy can be derived 
from the one on an appropriate dual loop algebra. 

Since our matrices A^ n > depend on z and A and the variable z appears both implicitly, through w{z) and 
its derivatives, and explicitly, we need to introduce some notation. Given any function / of A, z,w, w z , . . . , 
let us denote the partial derivative of / w.r.t. z as follows: 

* * 3/ a/ df 

OzJ ■= T7- + -^— W z + q w zz + ■ ■ ■ , 

oz aw ow z 

and dff the partial derivative of / considered as a differential polynomial of w depending on z, A: 

df 

QW f .— _L 

zJ dz 

Analogously d\f denotes the partial derivative of / w.r.t. A. Then given the matrices B and A^ as in (|I5|I 
and JIEJl, one has 

(21) d™A {n) = d x B, 
so that equation (|17|l is equivalent to 

(22) (d z -d™)A {n) = [B,A {n) ]. 

Remark 4.1. This phenomenon, i.e. equation \21\) . is a common feature of all Painleve equations and it 
was used in [2'6\ to find the algebro-geometric Darboux coordinates for the six Painleve equations. As far as 
we know, a proof of the fact that \21)) is a common feature of all the isomonodromic deformations equations 
is still missing. 

We are now going to interpret the evolution along (d z — d™) as a vector field on a coadjoint orbit of an 
element of an appropriate twisted loop algebra. Let LG be the group of smooth maps / from S 1 to SL2 
such that 

' 1 



/(A)a 1 (/(-A))- 1 =7, <r 1= ( 



1 



and A is considered as a parameter on S 1 . Denote by L271+2G the subgroup of maps of the form / 
/ + A~ 2 ™~ 2 /oo, where /oo is holomorphic outside S* 1 and let 02n+2 be its Lie algebra: 

{-2n-2 ~j 
X{X)= XiX^Xi e st(2,C), X(X)cr 1 = o-iX(-\) > . 

Then let G be the quotient of these 2 groups 

G = LG I L,2n+2G. 

Its Lie algebra is given by 

= I X(X) = J2X 1 \ 1 \X 1 e 0[(2,C), X(A)<7a = <tiX(-A) I / 02n+2 , 
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with Lie bracket defined as: 

[X{\),X{\)] = [ J2 [X k ,X i _ k ]\x i modg 2n+2 , 



oo / z+2n+l 



i=-2n-l \fe=-2n-l 



which obviously gives [X(A), X(A)] G g and satisfies the Jacobi identity. The dual space g* can be identified 
with 

9* = |H(A)=£H i A i |iVeN,S i e S t(2,C),H(A) < 7 1 =-(r 1 H(-A)|/^ n+2 , 
\X(X)= jr X i X i \X i €sl(2,C),X(X)a 1 = a 1 X{-X)\. 



02n+2 



2n+l 



by the following pairing 

(23) (X(A),H(A)) :=Tr(ResX(A)3(A)), VX(A) e g, H(A) £ g*, 

where Res indicates the formal residue, i.e. the coefficient of the A -1 term. Consider the subalgebra 

(24) g_ = S^X(X)=J2x i X i \X(X)a 1 =<T 1 X(-A)|/ fl2n+2 , 
its dual space can be identified with 

(25) g*_ = js(A) = f^A*^ G «[(2,C), S(A)a 1 = -<r x S(-A) | / fl2n+2 . 

An element X in the Lie algebra g acts on an element 5 £ g* by the coadjoint action 

(26) (ad^S, Y) := - (S, [X, Y]) = ([X, 3], Y) 
for any 7 eg. This shows that for every X G g, 5 6 g* 

[X,H] = ad^~ Gg*. 



When we restrict the coadjoint action to the subalgebra g_ and to its dual space gl, we obtain the following 
identification 

(27) [X_,H]+=ad^_H, S € fll, X_ e g_ 

where (•)+ is the projection from g* onto gl and (•)_ denotes the projection onto g_. 
Lemma 4.2. Given the matrices B and A^ as in 1^15)) and \lb)) . one has 

(28) [B,A in) ]=ad* B A, 
'A™ A" 2 " +1 



where B = 



Proof. Using 115(1 and 1(16(1 . we notice that: 

Then using the Drinfeld-Sokolov trick: 

[B,A {n) ] = 



G g_ and A = (A^*) + G Q*_, which is the dynamical part of A^ r ' 

_4(™)^-2"+l ' 



_4(«) \-2"+l 



,A (n 

+■ 

(^>) + 

where the last step is due to the fact that ^ — — — ^ commutes with (.4 



(29) 



/ >A (r i ) A -2„+l 
/ iA („) A - 2 „+l 

V 4™ 



□ 
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Remark 4.3. The fact that we can neglect the singular part of A^ at A = is more general than in the 
above proof. Suppose we want to compute a.d Xk A (n> for Xk = — — p ) , then for every Y G g_ we 



h 



ave 



• /^™)A- 2fe+1 "> 



2k+l 



4 A' 

because the singular part of A^ does not contribute to the residue. 
Similarly it is easy to see that 

{d z -d™)A {n) = (6,-8?) A. 
Resuming, we proved the following Lemma: 

Lemma 4.4. The monodromy preserving deformation equation jJ7| ) is the same as 

(30) (d z — d%) A = &d* B A, 

where A — (A^) + 6 0* is the dynamical part of A^ n \ and B — ^ 2 " j £ g_. 

This Lemma allows us to interpret the evolution along (d z — df) as a vector field on a coadjoint orbit of 
the twisted loop algebra g_. 

Let us now recall the Poisson structure on gl. This is fairly standard (see for example the beautiful book 
but we recall some details here in order to fix notations and adapt the computations to our special case. 

The Poisson structure on is given by observing that every X G g_ defines a linear function X* on 
9- 3 2: 

v . 0* ^ C 

This fact allows one to identify 0* with g_ and to define the Poisson bracket between two linear functions 
on q*_ as S. Lie did 

{X„n}(3):=(S,[X,Y]>. 
The Poisson bracket between two functions / and g on g* is given by 

(31) {/,<?} (H) = (3, [d/,d<7]) 

where the differential df of a function / on q*_ is a linear function d/ G ~ 0- defined by 

(32) (df, iH> := / (S + 6 X E) - f(E) + 0(6 X E) 2 , 
where 

5 X E := ad^S G 0* . 

In particular one has dX* = X. 

It is well known that this Poisson bracket is degenerate and its symplectic leaves are the coadjoint orbits 
of its elements. The kernel of this bracket consisting of the Casimirs, i.e. functions / such that 

ad^S(d/) = 0, VXg -,SG0*_. 

Lemma 4.5. The times ii, ... ar & the Casimirs of the Poisson bracket 

Proof. We show that the Hamiltonian vector fields generated by t\, . . . , t n -i are zero. 

Denote the eigenvalues of A(X) by ±/i(A). The polynomial part of /i(A) is a polynomial of order 2n and all 
the coefficients of the odd positive powers of A are zero. Therefore the polynomial part of /i(A) has exactly 
only n + 1 non zero coefficients, of which the first one is 4". Our claim is that all the other coefficients of 
the positive powers of A give our times t\, .. . ,t n -i (we shall see in Section [7| that the coefficient of the —2 
power of A is the Hamiltonian in the variable z). 

In fact, as proved in corollary 18. 21 below, the times U are given by the 'spectral residue formula' as 
follows. 

U = lRe Soo A- 2i - V(A)dA = {A, A-^vfr^-i, ) 
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where * is the eigenvector matrix of A and 03 = f J ^1 ) Since A 24i ^a^^ 1 commutes with -4(A), 
it produces a trivial vector field and therefore ti is a Casimir. □ 
Thanks to the above Lemma, the coadjoint orbits are obtained by fixing the values of t\, . . . , f n _i: 

{/ „(") \2fc V^n ?.(«) \2fc-l j_ \2fc N 

I l^k=0 a 2k+l A l^k=l°2k A + l^k=Q °2k+l A 

I v^™ J,(")\2fe-1 ri"-l \2fc „(") \2fc 

\2^fc=l°2fc A ~ 2^fc=0 °2fc+l A -2^fc=0 a 2fc+l A , 

The dimension of the coadjoint orbits is 2n. It is well known that the Poisson bracket l|31[) restricted to the 
coadjoint orbits is non-degenerate, so that it defines the so-called Kostant-Kirillov symplectic structure to on 
them. We will write 

u(f,g) = {f,gh 

for every pair of functions on the coadjoint orbit. 

To compute the Poisson brackets between the coefficients c4fc+i) b 2k ano - b 2k+ii we observe that their 
differentials are 



da^i = \{Eu -E 22 )\- {2k+1 ^ for 0<k<n, 



1 2fc+l 

(33) db^ +1 = ^(-E 12 + E 21 )\-<- 2k+ ^ for Q < k < n - 1, 

db 2k = l(E 12 +E 21 )\- 2k for 0<k<n, 

where with a slight abuse of notation we are calling a^fe+i' b 2k ano - b 2k+i ^ ne elements of jjl* ~ g_ which 

applied to E produce the coefficients ajjfc+i' &2*? an d ^2fe+i respectively. 

By using these gradients, we can compute the Poisson brackets between the matrix entries 

{ a 2 n k\v b ^li} = - h 2(k+i+iy for 0<k<n, 0</<n-l, k + l<n-l, 

(34) {a^k+v^} = - b 2(k+i)+i for < fc < 71, 1 < I < n, k + l<n-l, 

{ b 2 Mm} = 4^+0+1 for !< fc <^ 0</<n-l, k + l<n, 
while all the other brackets vanish. 

5. Canonical coordinates for the isomonodromic deformations 

Our first attempt to build the canonical coordinates for the second Painleve hierarchy is to use the general 
framework of the algebro-geometric Darboux coordinates (see ^H|j|SlE]). In this setting one considers the 
spectral curve 

(35) r(/j, A) = {det(/x - ^ (,l) (A)) = o} = {/i 2 = - dot (U (n) (A))} . 

The characteristic equation p? — — det (A^(X)) defines the eigenvalue /j(A) of ^4(A) as a function on the 
corresponding 2-sheeted Riemannian surface of genus g. The Baker- Akhiezer function "0(A) is defined then 
as the eigenvector of A^ (A) 

^W(a)v(a)=m(a)V(a) 

corresponding to the eigenvalue //(A). Generally, tp has g + 1 poles. 

Following |37| , let us briefly illustrate how to construct canonical coordinates p\ , . . . , p g and q\ , . . . , q g on 
the cotangent bundle of the Jacobian T* J of the curve T. 

Denote by q the A-projection of the generic point in the divisor of ip. We fix the following normalization 

(ci,c 2 ) ■ ip(q) = 1, 
for some choice of c\ , c 2 . The qj variables are the roots of 

(36) c\A 12 {q ) - Cl c 2 (An(qj) ~ A 22 { qj )) - c\A 21 { qj ) = 0, 
while the pj variables are the eigenvalues 

(37) p j =(A n (q j )--A 12 (q j : 

V C 2 
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Choosing the normalization ci = — c 2 = 1 we get roots q\, . . . , g 2Tl such that q n +j — —qj, j — 1, . . . ,n. They 
are the roots of the following equation: 



n-i 

(38) E^ti+4t^ 2k +^ 2n 

fc=0 

The corresponding pj are given by 

k=a 

In the generic case, it is well-known that the coordinates q\, . . . , g 2n , Pi, ■ ■ ■ ,P2n are canonical with respect 
to the Konstant-Kirillov Poisson structure as well 

(39) {pi , p 3 } = {qi , qj } = 0, {pi , qj } = % , 

however, a proof of this fact for non-generic cases is still missing. 

Generically, the dimension 2g of T* J coincides with the dimension of the symplectic leaves in the coadjoint 
orbit associated to (|30|l . This allows one to identify these symplectic leaves with T* J and to treat pi, . . . ,p g 
and qi, . . . , q g as canonical coordinates on the symplectic leaves themselves. 

In our case instead, it is not hard to realize that the suitably de-singularised spectral curve T is an 
hyperelliptic curve of genus g = 2n, so that dim(T* J) = An, which is twice the dimension of our symplectic 
leaves. In fact the characteristic equation has the following form 

/i 2 = 4 2 "A 4n + P / 2 „-i(A 2 ) + ^|, 

where Po( 2n _i is a polynomial of degree 2n — 1. By doing if necessary a small monodromy preserving 
deformation, we can assume this polynomial to be irreducible. Setting ft — Xfi, we get 

f? = Po/ 4rl+2 (A), 

where Pol^n+2 is an irreducible polynomial in A of degree An + 2. We see that the genus is 2n. 

Another problem is that the coordinates qi , . . . , q g are defined by taking the roots of the polynomial Ij38(l , 
so they may not satisfy the Painleve property of the isomonodromic deformations equations (see \A'2\ 
Therefore we propose a new set of canonical coordinates: 

Theorem 5.1. Consider the following 



(n) k( n ) n 
_ a 2(»-fc) + l + °2(n-fc) + l n V- 1 ,(n) 9S 2 j 
r fi — il 2fc — > l^fc — ~°2i ~arr ' K — i, 



a 2n+l j-- 



1 2j -■' mi,, : 



,,n, 



where Sk — Y)j—i Qj f or k = 1, . . . ,2n and Eli, . . . ,Il2 n are the symmetric functions ofqi, . . . ,q2n- 
ni=<7i+g 2 H hqin, n 2 = ^2 VjQk, n 2n = qiQ2 ■ ■ -?2n- 

l<j<k<2n 

Then 

(1) Pi, . . . , P„, Qi, . . . , Q n are coordinates in the symplectic leaves. 

(2) Pi, . . . , P„, Qi, . . . , Q n are canonical, namely 

{P, Pj} = {Q,, Qj} = 0, {PuQj} = 6ij. 

Remark 5.2. Observe that the first statement of our theorem could have been guessed by noticing that, as in 
the case of the Kowalevski top, T admits one extra symmetry A — » — A apart from the hyperelliptic involution 
n — > — fi. As a consequence T is a two-sheeted covering of a genus n hyperelliptic curve C obtained by setting 
z = \ 2 : 

C = | M 2 = 4 2 "z 2 " + Pol 2n ^(z) + ^ 

Differently from the case of the Kowalevski top, this cover is branched at and at oo. Having fixed the 
normalization c\ — — c 2 = 1, we see that the Baker- Akhiezer function ip has a simple pole at A = oo and 
2n simple poles at X = q±, . . . ,qi n . This pole divisor is clearly invariant w.r.t. the involution A — > —A. In 
fact the poles A = gi,... ,g 2n come in pairs qj,q n +j and each pair projects to one pole on C. Due to the 
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construction by A. Weil (see for example |H4j ). the symmetric functions of the pole divisor on C appear 
naturally as coordinates when endowing Jac(C) with the structure of algebraic variety. 

Proof. The first statement of the theorem follows as a straightforward corollary of Theorem 16.11 proved in 
Section El below. 

Let us prove the second statement of our theorem, i.e. that our coordinates Pi, . . . , P n , Q\, . . . , Q n are 
canonical. 

To compute the Poisson brackets between our coordinates, we observe that H 2 k+i = S 2 k+i = in our 
case and 



» 



n 



2 k 



a 2{n-k) + l + ^2(ri-/c) + l 



An) 
l 2n+l 



First let us compute the bracket {Pk, Pi} 



+ h yn> n (n) + b ( 

a 2(n-fc) + l T U 2(n-fc) + l u 2(n-/) + l T u 2(n~l) + l 



k = 1, . . . , n. 



dn) 



{Pk,Pl} = 



(n) 
I 2n+1 



An) 
l 2n+l 



u 2n+l / 

1 \ 2 
u 2n+l / 



,(«) 



» 



°2(n-l)+V fl 2(ri-fc) 



dn) 



»}-{ 



,(«) 



An) 



^2(n-fc)+l' a 2(n-l) + 



.}) 



2(2n-2-fc+l) w 2{2n-l-k+l) 



Therefore we have 
(40) 



{P k ,Pi} = {U 2k ,U 2 i} = 0, k,l = l,. 



To compute the brackets that involve the Qk, we make use of the following formula: 

N fe-l»5fc _ k 



(41) 



fc=i 



fc" 7 ' 



where 7 is an auxiliary variable. Since in our case, the roots of the polynomial 



E( 6 2r + i+4ti)A 2fe +4v + iA 2 " = o 



fc=0 

are given by qi, . . . , q n , -qi, . . . , -q n , we have 

Il2fe+1 = &2k+l = 0. 

Therefore, by differentiating il4*T|) with respect to n 2 fc, we can express ^j^- as follows 

dS 2j 



(42) 



an 



2 A- 



"(2j) 



\i=0 



2j-2k 



where [JT(7)]2j_2fc is the coefficient of -f 2 i~ 2k of ^(7) considered as a power series in 7 near 0. Note that 
although the sum in the left hand side of (|41|l goes from 1 to 00, only terms where j < n enter in (|42(l as 
j — k < n. 

We will now compute the bracket {P k , Qi} 



(43) 



{iVQi} = £{n afc ,&$} 



(«)-, 1 dS 2 j 



2j 8U 2 



2,/ 



2j 



an 



2/ 



We are grateful to M. Talon for pointing out to us the similarity with the Kowalevski case and with the construction by 



A. Weil. 
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Since gjp- is a polynomial in n 2m with m < j, the second term in l)43|) is zero because of l|4Uf) . The first 
term in 1)43(1 is 



2j an 2 . 



E n 2C=-j) 



E n 2(fc-j) 
3=i 



E n 2/ 



2.; 



2^-2; 



2j-2l 



where we replaced the sum from 1 to n by a sum from I to n in the last equation because the expression 



does not contain any negative power. Since the expression 



E n 2(fe-J') 
3=i 



E n ^ 2j 



2j-2i 



is just the coefficient of r y 2k 21 in 



A 2iT" I I E 11 ^ } = h 
■,'■■<< ) \j=0 



E n 2 



2i 



we see that 

{P k ,Qi} = S kl . 

To compute the brackets between Qk and Qi, once again, note that since the brackets 

1" 2J A, >-\an 2fc 'an 2 J" u ' 

the only contributions to the bracket {Qk, Qi} come from the cross terms 



{Qk,Qi} 



E 



2j 2j an 2fc ' 2t ' 



2i 



an 



(44) 



j, i=i J 



(n) 
2j 



95; 



an 



2i h (n) 
J °2i 

2k 



2^ 4 y I °2, 



as- 



an 



2j_ .(n) 
' °2i 



2 A' 



2 j 1 








hi J 








552, 

an 2Z 


^< 


fas^ („)' 
ian 2/ ' 2J J 


as 2j 
an 2fc 


dS 2l 


4?] 


(„)' 

ian 2i ' 21 J 


I os 2i 


an 2 . 


on 2k 
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The bracket between q^ 3 and bi^ can be computed 



as 



<• ZK ' \m=0 / s=0 

2 



7 



2 s 



= -2.7 



-2.7 



= -2.7 



= -2.7 



2j-2k 



2(s-i)7 



\m=0 



s=0 



EwH E n 



-2s 7 



,2s+2i 



\m=0 



s=0 



2j-2fe 



2j-2k 



£n 2ro7 H E n 



2s 7 



\m=0 



s=0 



EW m E n 



2s7 



\m=0 



s=0 



2j-2fc-2i 



i- -i 2j-2fc-2i 

where in the last line, we replaced the upper limit of the second sum by n. This is because 2j — 2k — 2i < 
2{n —i), so if p > n — i, the coefficient of j 2p in 53s=o n-2s7 2s will not enter in the final expression. 



Therefore we have 



1 fftStf An)} ds 2i _ 



4ij \ du 2k 



' u 2i 



dn 



21 



-3 



EW m E n 



-2s 7 



2.9 



\m=0 



s=0 



Similarly, the second term in l|44|) is given by 

El ( dS 2j dS 2i 

therefore the first and second terms in (|44ll equal each other and 
In summary, we have 



2j-2fc-2i 



EwH E n 



-2s 7 



Vm=0 



s=0 



2j-2fc-2i 

we have {Q k , Qi} = 0. 



as we wanted to prove. 



{P k ,Qi} = S u , {P k ,Pi} = 0, {Q k ,Qi} = 0, 



□ 



Example 5.3. Case n = 1. In this case Q\ = ^^fn 2 " = — an< ^ Pi = II2 = Ql pj 1 , that is 

2 a 3 



These 



Q = 4w, P - - (to, - w 2 - ~ 
coincide with Okamoto's canonical coordinates (up to a constant factor) . 

a< 2) +6< 2 ' 



Example 5.4. Case n = 2. In this case Pi = II 

'2 



°3 + °3 p _ TT 
(2) 1 r 2 — 11 

i & l 2) |f| - -4 2) + ^ 2) n 2 , and Q 2 = i&ffgj = -bf so, finally 



(2) ! ^1 — 2 " 2 gn 2 



p 



1 



2 — \-z + 6w 4 - 12w 2 w z + 2w 2 z - Aww zz + 2w zzz + 2t 1 (w z - u> 2 )), 
lb 

Qi = -8ww z + 4w zz , 



Q 2 = 16w. 

6. The coefficients of A^ as polynomials in the canonical coordinates 

This Section is completely devoted to the proof of the theorem below which expresses the matrix A^ n ' as 
a polynomial in P x , . . . , P n , Q 1 , . . . , Q n and z, ti, . . . , i n _i. 
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Theorem 6.1. Let A, B oc id and B even be the following polynomials in X: 

n n n 

(45) A = Y, a %l X2l+1 > ^ = E^iiA 2i+ \ B even = Y,b { 2 l } \ 2 \ 

i=0 i=Q i=l 

and V , Q and T the following truncated series: 

n n n — 1 

(46) Q = ^Q,A 21 , V = Y,P*^ 2 \ T = 5>(2A) 2l - 2n -z(2A)- 

i— 1 i— 1 i—1 

Then the following relations hold true: 

1 /o\>2n+l (-> , tj _ ( 1 + T f\ _ loY ,-2n-l n 2t 



(47a) A = l^-(2\y n+L -(2\)- Zn - L Q z (l + T) 

(47b) B odd = Q(2A) 2 " +1 U+V- ^|r ) +(2A)- 2 "- 1 Q 2 (l + 7^) 

(47c) £ eue „ = -A 2 (A- 2 Q(l + 7')) + , 

where in the above, the inverse (1 + is £o 6e interpreted as: 

oo 

(48) (i + ^)-i=^(_p)i 

t=0 

and X + indicates the polynomial part of the Laurent series X . 

In particular, by comparing the coefficients of A in j47\) , we can express the matrix entries a^^, ^i+i 
and 62™' as polynomials in the canonical coordinates and the times. 

Proof. Throughout the proof of this theorem we shall use the following facts: let X, X' and Y, Y' be Laurent 
series in A with no positive part, such that 

(\ 2n ~ 1 X) + = (A 2 "- 1 X')+ 

then 

(49) (A 2 "- 1 iy)+ = (\ 2n - 1 X'Y')+, and (A 2 "- 1 ^ 1 ^ = (A 2 "- 1 (A:')" 1 )+. 
Let us prove (|47c|) first. Using the definition of Q k given in Theorem 15. II we have 

An) _ -I ( n ) dS 2k _ 



2k 2k 2k dn 



2k 

1 An) 9S 2 k+2 1 An) dS2k+i 

2k + 2 2k + 2 dn 2k + 2fc + 4 2k + 4 dU 2k 



(n) 

now apply the above repeatedly to b 2k+2 i ^° OD tain the following 



6g? = -Q k -Q k+ J 1 U ° 2k+2 



(50) - Q fc+2 



where 



1 dS : 
2kT2~dh 2k 
1 dS2 k +4 . 1 1 dS 2k +4 dS 2k+2 



2k + 4 dU 2k 2k + 4 2k + 2dn 2k+2 dll 2k 
Qk+l ( oh I 9.1 5 fnT' + U ' 2k + 21 ) + " ' 



m— 1 



/-^ r/m = \ -* 1 dS 2k+2 i \ - 1 dS 2k+2jl 

1 j 2fc+2 ' A 2fc + 2; an 2fe+2jl . ^ 2 fc + 2* dn 2k+2n 



V - 



dS 2k+2jm _ 1 1 dS 2k+2jn 
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Note that, by (02), we have 



(52) 



dS- 



2p 



on 



2q 



"(2p) 



vJ=0 



2p-2q 



This allows us to interpret the t^^ +2 i as coefficients of an infinite series, which is a product of the series in 
the right hand side of 1)52(1. However, we need to be careful since the sums in (|51H begin with terms that 
correspond to the coefficient of 7 2 in (|52[l . This means that 

m+1 



/ 



^2fc+2Z — coefficient of 7 2i in 



1 



V 



E n 2r 



2.; 



We can now identify Qk as coefficients of A in the polynomial Q, b 2k as coefficients of A in the polynomial 
B even . If we now substitute 7 = A -1 and interpret U%l +2 i also as coefficients of A in the product series, we 
can rewrite 11511) as follows 




which proves l|47c|) . 

Let us now prove l|47b() . 



Notice that in the expressions (|15c|) for b 2k+1 the Lenard operators appear 



together with their first and second derivatives. The canonical coordinates instead involve only the Lenard 

(n) (n) 

operators and their first derivatives. In fact, the P±, . . . , P n are expressed in terms of sums a^j+i + ^2j+i 
which depend only on the Lenard operators without derivatives, and the Q\, . . . , Q n depend on the same 
sums and on the even coefficients b 2k which are only expressed in terms of the Lenard operators and their 
first derivatives. 

Therefore, to express the odd coefficients 6 2 fe+i m terms of the canonical coordinates, we will need to 
express the second derivatives of the Lenard operators in terms of the Lenard operators themselves and their 
first derivatives. In terms of the generating function, this is given by 

Proposition 6.2. Let C be the generating junction of the Lenard recursion operator: 

00 

C := J2 Ci ?> 

i=i 

where £ is an auxiliary variable. Then the following relation holds true: 



(53) 



1 



#C = r(r^-A)-2£A 



(l + 2£)- 



In particular, the above equation expresses the second derivatives of the Lenard operators as a polynomial in 
the Lenard operators and their first derivatives. 

Proof. We prove (|53|l at each order in £. At order £, the equation (|53|l is trivially satisfied. At order £™, 
n > 0, by multiplying Ip^)) by (1 + 2£), we get: 

(54) C n+1 = dlC n + iC n C l + ^(^C l C 3 C n -, J -d z C, J d z C n - J -C 3+l ^^ 
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To prove this, we integrate relation by parts iteratively. At the first step, we get: 
(55) £„+i = (d 2 + 2£ x ) £ n + 2 J dd z £ n dz, 

where we have replaced w z — w 2 by £\. It is therefore sufficient to compute J £\d z £ n dz. To achieve this, 
we will now compute a more general term J Cid z Ckdz for any k and i. First let us replace the term d z £ k 
by using (|T3|l : 



J dd z £ k dz = J d (d^ + A£ 1 d z + 2(£ 1 ) z )£ k _ 1 dz 



We now integrate the first and last terms by parts, where the integration by parts in the last term is performed 
as follows, let 

F(z) = 2 J d{£ x ) z dz 

= dlLi + ALxd-d+x, 
where we used the Lenard recursion relation Ijl^J) to obtain the above equality. We then have: 

2 J A(£i)*£fc-idz - F(z)£ fe _i - J F(z)d z C k - 1 dz. 

; From this we have 

dd z £ k dz = dd^Ck-x - d z Cid z C k -x + £ k -id 2 d + 4£iA£fc-i 
— d + i£ k -i + J £i + \d z £ k -idz 

We can replace the term J d + id z £ k -\dz on the right hand side by similar formula and express J dd z £ k dz 
as a polynomial in £j, d z £j and d 2 £j for j < k. In particular, we have 

2 J £id z £ n dz = d£ n + ^ ^4£ 1 £ 3 £ n - j - d z £jd z £ n ~j - 
£j+i£n—j 2£jQ z £ 7 i— j \ . 

By substituting this into 1)55(1 . we finally get (|54|l . □ 

in) 

We can now express the odd coefficients b 2k+1 in terms of the canonical coordinates. 

If we make the substitution £ = \X~ 2 in the generating function £, then by l|15c[> . we have 

(56) B odd = ((2X) 2n - 1 (ld 2 £ + wd z £ + w z (£ + £ ))(l + T)] , 



2 

(57) B even = (2X) 2 ((2X) 2n - 2 (-d z £-2w(£ + £ ))(l + T)) + . 

(58) 

Note that, since 

((2A) 2 "- 1 ^£(1+T)) + 
is a series in odd powers of A only, while 

((2X) 2n - 2 d z £(l+T)) + 
is a series in even powers of A only, we have 

2A {{2X) 2n ~ 2 d z £{l + T)) + = ((2A) 2 "- 1 ^£(1 + T)) + . 
^From this, and the second equation in (|56|l . we can express the second term in B od d in terms of B e 
((2A) 2 "- 1 a 2 £(l + T)) + = {2Xy 1 B evm + ((2A) 2 ™- 1 (2£ + l)u,(l + T)) + . 
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By substituting this into the right hand side of B odd , we have 

B odd = ^X) 2 ^ 1 (^d 2 z C-w 2 (£ + C ) + C 1 (C + Co^j(l + r)^(2X)- 1 B ev 
B even = (2A) 2 ((2\) 2n - 2 (-d z £.-2w£-w+)(l + T)) + . 

By substituting this into (|53J) . we obtain 



B odd = ^((2A) 2 



2w 2 L-w 2 + i((2A) 2 £ + A) 



(59) 



(1 + 2Cy 1 ((2A) 2 £ - A) - CCi 



(1 + T) -2w(2X)- 1 B e 



To complete the calculation, we need to express C in terms of V . Since we have 



» 



Pk = 



2(n-k)+l T "2(n-fe)+l _ 2 ^ 

~ 4fe 2^ Ht-k+V-n), 



An) 
l 2n+l 



k = 1, . . . , n, 



l—n — k 



where i„ = 1, we see that 



£fc — 2 2k l P k — t„_i£fe-i 



C k = 2 2k - 1 P k + 2 2k - 3 P k - 1 (-t n - 1 ) + 

(n—i—l 
-<n-i+-+ E K-HJEK-i,)- E (-*«)• 
ii=l 32=1 j, = l 



n-k-1 

(~^n— fc— ji ) ^ 1 

Jl=l 32=1 J(=l 



*)••■ E Hi,) 



Finally, by using similar argument as before, we see that 

(60) ((2A) 2 - 1 (£ + Co)) + = ((2A) 2 "" 1 (l + ^) (1 + T)" 1 ^ , 

where (1 + T) _1 in the above is to be interpreted as in 148|) . Thanks to (|49|l . this implies the following 

(61) ((2A) 2 - 1 (£ + £ )- 1 ) + = (^A) 2 - 1 (a, + ' (1 + T)) 

By substituting (JgOJ) and ((HTJ into we get lj47b|l . 

To finish the proof of the theorem, note that by the definition of the coordinates Pi, . . . , P n in Theorem 
15.11 we have 

X(2\) 2n (l + P) = (A + B odd ) 



(62) 

The proves the theorem. 



A=(-B odd + X(2Xy n (l + V)) 



□ 
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Example 6.3. Let us illustrate how formulae (|47a|l . 147bjl and (|47c|l work in a concrete but non trivial 
example. In the case n = 3 the definitions (|45|l and (|46l) read: 

^4 3) A + 4 3) A 3 + 4 3) A 5 + 4 3) A 7 , 

B odd = b?h + bfh* + b?h 5 , 

(63) Q = QiA 2 + Q 2 A 4 + Q 3 A 6 , 

P= P +^+^ 
A 2 A 4 A 6 ' 

j- h , ti z 



(2A) 2 (2A) 4 (2A) 6 ' 

Theorem 16.11 allows to express the coefficients a^\a^\ dj?\ c^p and the coefficients . . . , fog 3 -* in terms of 
Pi, P 2 , P3, Qi, Q2, Q3 and z, ti,t 2 - Let us use (|47c|) first: 

Q(l + V)= + P2Ql + P3Q2 + PlQl + p 2 q 2 + P3Q3 + A 2 (g 1 + P X Q 2 + P 2 Q 3 ) + A 4 (Q 2 + P X Q 3 ) + A 6 Q 3 , 

and after dividing by A 2 and throwing away all negative powers, we get 

(A- 2 Q(1 + V)) + =Q X + PiQ 2 + P2Q3 + A 2 (Q 2 + P1Q3) + A 4 Q 3 . 
Now, we need to multiply by A 2 again and to compare with B even . We get 

4 3) = ~Qi - P1Q2 - P2Q3, bf = -Q2-P1Q3, b^ = -Q 3 . 

Let us briefly illustrate how to obtain the other coefficients. The procedure is the same as above with the 
only complication of the term 

Let us compute this term explicitly. Since we have a power 7 in A in the numerator, we need to only the 
terms up to order in the expansion of (1 + P) -1 at 00: 



(i + vr 1 = 1 - % + + 2PlP2 ~/ 3 = P! + o<n 



In this way we see that 

I(2A) 7 -j-±|^ = (64Pi - 16i 2 )A 5 + (64P 2 - 32P 2 - 4tj + 16t 2 Pi - 2i 2 )A 3 + 

+A(64P 3 + z + 32P 3 - 64P!P 2 + 4t ± I\ - 16i 2 P 2 + 16< 2 P 2 - ht 2 + 2f|Pi). 
Analogously, one can compute the other terms in (|47a|) and (j47b(l to obtain: 

6( 3) = z + 32P 1 3 -64P 1 P 2 + 64P 3 + ^i ' 0l ° 3 ■ PlQ2 ° 3 ■ 



128 64 64 

+ 1§ + 4tlPl ~ 16hP ? + 16t2p2 ~ hh + 2 * 2p ' 

b& = 64P 2 - 32P 2 + 9j9l + ^ _ 4ti + 16i 2 P! - 2i 2 , 
6f = 64Pi + ^| - 16* 2 , 

„(3) _ „ , 09 p3 aA p p , Ql , QlQ3 P1Q2Q3 

+ P rSr + 4t i p i ~ 16t 2^i 2 + 16i 2 P 2 - tit 2 + 2t|Pi, 



128 



4 3 > = -32P 2 + ^ + ^ - 4 tl + m 2Pl - 2tj, 

4 3) = +|| - 16t 2 , 4 3) - 64. 

It is clear that we have only used linear algebra to obtain these coefficients, instead of using differentiation, 
integration and recursion. 
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7. Hamiltonians 
In Section 0] we proved that the equation (|16|) is the same as (|3(J[1 : 

{d z -d™)A = Bd* B A, 

where A = (A^) + e g*_ is the dynamical part of A (n \ and B = (^-^r—^ G 0-- This allows us to 

interpret the evolution along (d z — Of) as a vector field on a coadjoint orbit of the twisted loop algebra g_. 
We are now going to show that this vector field is Hamiltonian and that the isomonodromic deformation 
Hamiltonian for the n-th equation in the PII hierarchy is given by 

(64) := -^TrRes ^A 1_2n (U (n) ) 2 j . 

This fact is actually a consequence of a more general result: 
Proposition 7.1. The vector field 



X k (A) 



is Hamiltonian with the Hamiltonian function 



A 



(65) 4 ,l) :=iTrRes(V- 2fe (.4( n ^ 2 



Proof. Let us denote 

L k = - (a^X 1 

We are interested in the vector field 



X k (A) 



L k ,A 



To show that it is Hamiltonian and to compute the Hamiltonian function /, we use the following definition: 

u{X k ,Y)(E) :=-([Y,H],d/), Y e fl_, 3 S fl* , 

so that 

u(X k ,Y)(A) =-([L k ,Y] ,A) = {[Y,A],-L k ) =< [Y,A],df > . 

This shows that if we can prove that there exist / such that d/ = —L k , then [Lfc,A] defines a Hamiltonian 
vector field of Hamiltonian /. 

We are now going to show that the Hamiltonian (|65|) is such that d/u — —L k . For every l£g_ and 
S € Jjl, we can identify [X, S] = ad^S with a vector tangent to the coadjoint orbit. Denote this vector by 
S X S. Let SxA £ TaOa then using the definition 1321 we get 

h£\A + 5 x A)-h k n \A) + 0{8 x Af = l -(\ l - 2k A, (26 X A^)) = (s x A,dh k n) 

which is the contraction between 5 X A and dh k , as we wanted to prove. □ 

We now compute the Hamiltonian lHS n ) in terms of the canonical coordinates. 

Theorem 7.2. Define 
(66) 

/n— 1 n — 1 n \ ~ 

n ^!,...,x- n , m,...,\4n,Z) — 4n I 2(n-l)-l / , "2i+l"2(n-Q-] "*" 2f 2(n-Q j 4" ' 

\i=0 i=0 i=0 / 

in which we are thinking of c4™4-i) ^2™+i' ^2? as ^ e functions of Pi, ... , P n , Qi, . . . , Q n , ti, . . . , z given 
by [47j )- Then the n-th member of the second Painleve hierarchy is given by the equations 

( 7) dPk = &HW dQ k = 

1 ' dz dQ k ' dz dP k ' 
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Proof. Observe that equation 1(64(1 gives 

^ / n—l n — 1 n \ 

- ff< " ) = _ 4^7 ( X/ a 2J+l a 2(n-i)-l ~~ X/ ^2i+1^2(n-;)-l X] ^2Z ^2(n-i) ) ' 
\;=0 1=0 J=0 / 

where we are treating a^, 6^ as coordinates on the coadjoint orbit. When expressing this Hamiltonian in 
our canonical coordinates, we need to take into account a shift h due to the explicit z dependence in the 
variable P n . All other canonical coordinates depend on , bj only. To compute this shift we use the 
following well-known result: 

Lemma 7.3. Let 

(68) ^ = { Vi ,H(y,z)}, 

az 

be a Hamiltonian system on a Poisson manifold with Poisson brackets {•, •} and 

be a local diffeomorphism depending explicitly on z. Let the vector field d z <p be a Hamiltonian vector field 
with Hamiltonian SH. Then $68\) is a Hamiltonian system also in the x -coordinates 

^ = {xi,H(x,z)}, 

where 

H(x, z) = H(<fi(x, z), z) — SH(<p(x 1 z), z). 

Let us compute this shift in our case. The only coordinate depending explicitly on z is P n — + 
/(a£ n) , a { 2 n J +1 ,b { n \b[ n \ . . . , So for y = (P 1: . . . ,P n , Q 1 , . . . , Q n ), we have 6HW = %t which gives 

TO- D 

Remark 7.4. It is clear that the Hamiltonian equations j 6'7| J satisfy the Painleve property. In fact ji7| ) 
satisfies the Painleve property 32 33 , and since 

A (n) 

is a polynomial in Pi, . . . , P n , Qi, . . . , Q n , no movable 

critical points are introduced. 

Example 7.5. In the case n = 1, we get the Pn Hamiltonian (361 129j 

H (l) = 4 p2 + 1 g + lpg2 + 2 p z _ Ig^ 



The Hamilton's equations 



az* 

give the second Painleve equation for w = jQ: 

d 2 w 



d " Q -\Q 3 + Qz + A ai , 



, „ = 2w 3 + wz + a\. 
az z 



8. Higher order flows as time-dependent Hamiltonian systems 

As illustrated in Section |2 the time t k dependence is described by the rescaled mKdV equation i|14|) , 
which is equivalent to the compatibility equation 1(18(1 . In the Appendix we proved that all equations 
1)17(1 . 1)18(1 and 119(1 are consistent, so that they indeed define isomonodromic deformations. In this Section, 
we shall deduce the Hamiltonian functions W]l such that equation ((14|) is equivalent to 

1 ' dt k ~ dPi ' dt k - oqr i - i '---' n - 

Since we are dealing with isomonodromic deformations, the correct way to deduce these Hamiltonians is to 
express equation ((19(1 as a time-dependent Hamiltonian system. 
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8.1. Spectral invariants. The spectral curve T of the matrix A^ defined by 

detGitl- A (n) ) = 

plays an important role in defining the Hamiltonians of the isomonodromic flows and of the Painleve hierarchy. 
We now express T in terms of the generating functions of the canonical coordinates and the time: 

M 2 = \- 2 (A 2 -B 2 odd +(B even -a n ) 2 ) 

= ((2A) 4 "(1 + V) 2 - 4(2A) 2 "- 1 (1 + T)B odd + \- 2 {B even - a n ) 2 ) 

By substituting (|47|l into the above, we obtain the spectral curve in terms of the canonical coordinates 

H 2 = (2A) 4 "(l+7>) 2 -4(2X) 2n - 1 (l+P)B odd + \- 2 (B even -a n ) 2 = 

= (2A) 4 "(1+P) 2 -(2A) 2 "- 1 (1+P) 



4(2A)- 2n " i Q 2 (l + 7') 



(2A) 2 "- i (l + P) 



(2A) 



(2A) 2n +1 (^ l + p _(±±lpj 

A" 2 (A 2 [A- 2 Q(1 + 7')] + + c 

4(2A)- 2 "- 1 Q 2 (l + 7') 



l + V 



+ a- 2 (a 2 [a- 2 q(i + t')]_ 



where ^pp = Y^Lg(~^) 1 '^' 1 ■ The last equality follows be A 2n (l + V) contains positive powers only. 
In particular, we have proved the following 

Proposition 8.1. The coefficients of the spectral curve /i 2 of the matrix A^ can be expressed as polynomials 
of the canonical coordinates and the times as follows 



fi 2 = (2X) 2n - 1 {l+V) 



(2A) 



2n+ 



l( l+3~) 2 

l + V 



-A{2\)- Zn -'Q^{l + V) 



-i + 



+ A- 2 (A 2 [A- 2 Q(l + 7>)] -a n 



(70) 

Corollary 8.2. The constant a n , the times t\, . . . i n _i and the Hamiltonian are spectral invariants. 

In particular the spectral curve can be written as 



M 2 = (2A) 2 " ((2A) 2 "(1 + T) 2 ) + - 4 "ij(")A 2 "- 2 + £ h{ n) X 2k ' 2 + 



k=l 



A 2 



where 



» 



1 



-ResTr A 



,1-2/c 



A {n - 



k = 1, . . . , n. 



Proof. By definition of it is clear that h£ is the coefficient of the 2k — 2 power in A. Analogously 
for H^ n \ So we only need to prove that the coefficients of the powers 4n, 4n — 2, . . . ,2n are given by 
(2A) 2n ((2A) 2 ™(1 + T) 2 ) + . In particular, this implies that the polynomial part of fi is given by the following 

H = A n t n \ 2n + 4™- 1 t n _ 1 A 2 ™- 2 + ...-z. 
In fact, let the expansion of the spectrum /j, at A = oo be the following 

2n 



Let us denote the coefficients of fj, 2 by Di, that is 



i=-2 
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We can express Di as quadratic polynomials in the fi k and if k < 0, then ji}. will only appear in the coefhcient 
Di when i < 2n. Therefore, to compute the polynomial part of fi, we only need to consider the coefficients 
Di with i > 2n — 1. These coefficients are given by the coefficients of 

((2A)- 2 V)+ = ((2A) 2 " ((2A)- 4 V)) + 
By using the relations l|49l) in J7UJ, we see that 



((2A)- 2 > 2 )+ = 



(2A) 2n (l + T) 2 -4(2A) 



-2n-2^i2 



-I + 



+ [4(2A)- 2 "- 2 Q 2 (l+7>) 2 ] + 
((2A)- 2 V)+ - [(2A) 2 "(1 + T) 2 ] + 
This implies the corollary. 



□ 



8.2. Time flows Hamiltonians. We want to adapt the construction of Section0]to express equation II 19( 1 
as a time-dependent Hamiltonian system and the computations of Section to find the Hamiltonians. 
The main difficulty we encounter is that now 

d x M {k) (2k + l)d™A {n) . 

The main idea to handle this problem is the following: suppose there exists a set of coordinates ui, U2, ■ ■ . , U2n 
in our coadjoint orbit O a such that 



(71) 



d x M w = (2k + 1)9" A {n) 



where <9" fc denotes the ^-derivative with the u coordinates fixed (in the sense explained at the beginning of 
Section 2J. Then (|19|) becomes 



(72) 



(2k + 1) (d tk - dl) A (n) = [MW,AW] 



after cancelation. This allows us to interpret the evolution along d tk — <9" fc as a Hamiltonian vector field on 
the coadjoint orbit, as explained in Section 0] provided that the coordinates ui, . . . , U2n exist (which is a 
non-trivial fact because equation l|71|l is overdetermined). 

Our strategy is as follows: since in Proposition 17. II we computed the Hamiltonians corresponding to 
the matrices which give the same flow as 

(73) L k =\\ 1 - 2k A {n) ] , fe = l,...,n, 

L J + 

we introduce some new times s%, . . . , s n corresponding to the flows along L\, . . . .L n : 

(74) d Sk A^ = [L k ,A^]+d x L k , k = l,...,n. 
•^r. Since the matrices L k are related to the matrices M^ k > by 

n 

k = 1, . . . , n. 



L k = A k J2 tiM {i - k) 
i=k+l 



4 k t k B, 



Observe that s n 
(75) 

the times ii, . . . , t n -i and the si, . . . , s n -i must be related by: 

n 

(76) 8 Sk = 4* ( 2 ( z - k ) + i)^*.-. - 4 fe i fe 9 z , k = l,...,n. 

i=k+l 

This relation is the key to our procedure: on one side, it will allow us to prove that there exist coordinates 
ui j • • • > U2n such that 



(77) 



cT *4 (n) = d x L k , 



on the other side relation (|75|l , will allow us to compute the Hamiltonians H^' in terms of the Hamiltonians 

4 n) - 

The following proposition shows that l|7t)|l guarantees the compatibility of the over-determined system 

CU>- 
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Proposition 8.3. Let s k be the times corresponding to the equations 

d Sk A™ -d x L k = [L k , A™], k = l,...,n 
where L k are given by Then the system of differential equations 

8 u Sk A^ = d x L k 
is compatible only if t k and s k satisfy the following relations 

(78) d Sk t 3 = A k (2j + l)t j+k , d Sk z = -4 k t k , 
where ti — for I > n. 

Observe that equation l|78|) is equivalent to l|76[) . 
Proof. Let us consider the explicit s k derivative given by 

&: k A^ = o X L k . 

This is an over-determined system of equations. We are going to prove that it admits a solution. In terms 
of the entries of AS n > , this is equivalent to 

9s A% = (2j + l)ogJ. ab+1 , j = 0,...,n, 
d l h f ] = 3bj+M, i = l,..-,2n, 

af ] = 0, j > In + 1, bf ] = 0, j> 2n. 

This means that 

x 2n d: k (i + v) = [d x {\ 2n ~ 2k+1 {i + v))} + , 

(79) \~ 1 ds k B even = [d x (\- 2k B even )] + , 

y^Bodd = [dx(\~ 2k B odd )] + , 
where V, B even and B dd are given by ffify and l|47|l. We will now rewrite the expression of B dd in <|T7j) as 



Bodd — j 



(>>\\2n+lC\ , V \ ((2A) 2 " + 1 (1+T)) 2 



(80) + (2A)Sl(l+-p) 

Then, by computing 

(81) \- x d u Sh B odd = [d x {\- 2k B odd )] + , 

and applying lf79*|l and ijSUIl . we see that, in order that the over-determined system lf77)l is compatible, the 
derivative of T must satisfy 

X 2n dl{l + T) = [d A (A 2 "~ 2fc+1 (l + T))] + . 

This gives 

&l tj = 4 fe (2j + l)t j+k , 81 z = -4 k t k , 

where t\ = if I > n in the above equation. 
If we chose another set of coordinates, 

{yi(u,t), . . .,y 2n (u,t),ti, . . .,t„} 

Then the vector field t?f is given by 

n 

,r ;>■< . 

Therefore tj = d^ k tj for j = 1, . . . ,n and the same for z. Hence we can drop the superscript u in (|82() to 
obtain JTSJ). □ 
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As proved in Proposition the Hamiltonians /i£ n ' corresponding to the vector field 

X k ■= [L k ,A], 

are given by equation l|65|l . Due to our derivation, they must be thought of as written in the coordinates 
u. Since the canonical coordinates Pi, ... , Q n depends explicitly on the when expressed in terms of the 
coordinates u, the actual Hamiltonian differs by a shift given by Lemma 17.31 

We first compute the explicit Sk derivatives of the canonical coordinates. From the definition of the 
coordinates Q±, . . . , Q n we have 

X- 2 Q = -[X- 2 B even (l + V)- 1 ] + . 
We can the apply (|79|l to compute the explicit Sj derivatives of Q. This gives 

dl Q= [\ 1 - 2k d x Q-(2n+l)\- 2k Q] + . 
Hence the explicit Sj derivatives of the canonical coordinates are given by 

dlQj = (2(j + k)-2n-l)Q j+k , 

(82) dlPj = (2(n-j) + l)P i _* l jjkk 

dlP k = (2(n-*) + l), 

where if I > n or I < 1, then Pi — 0, Qi = 0. 

Using Lemma \l . 31 and equation l|82|) . we can compute the shifts : 

n—k 

(83) Sh{ n) = X)( 2 C? +k)-2n- l)PjQ j+ k + (2k - 2n - l)Q k . 

3 = 1 

Finally, we use the formulae above to show that the coordinates u actually exist. In fact we construct 
them recursively as follows using l|82() . For j = 1 this gives 

d^Px = (2n - l)S kl , 

so that we choose Ui = Pi — (2n — X)s\. Then for j = 2 equation l|82|) gives: 

8^P 2 = (2(n-2) + l)P!, 
&IP 2 = (2n-l) 

and all other derivatives of P 2 are zero, so we put 

s 2 

u 2 = P 2 - (2n - l)s 2 - (2(n - 2) + l)P x Si + (2(n - 2) + l)(2n - l)y, 
so that Ps k u 2 = 0. By repeating this procedure iteratively we get: 

U* = Pi+^WijPi-j +Wa, i = l,...,n, 
i=i 
n—i 

Ui+n = Qi + s.VjjQi+j, i=l,...,n, 
i=i 

where and Vij are given by 

W * = EE^-^ + J'-M-l)" 1 (f-) 1 ■■■(-2(n-i + j-l c )-l) k < 

771—1 \ / \ C / 

m=l \ 1 / \ c 

where the second summations in the above are taken over all possible combinations of integers 

hki + l 2 k 2 ^ \-l c k c =j 

ki + h k c = m, < ki < m, < U < j. 

Resuming, we proved the following: 
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Theorem 8.4. The Hamiltonians for the equations 

(84) d Sk A^-d x L k = [L k ,A^} 
where L k are defined by are 

(85) h{ n) = ^ResTr (V- 2fe (U^)^ 

and their shifts Sh^ corresponding to the canonical coordinates Pk, Qk defined in Theorem \5.1\ are given by 

n — k 

8h^ ] = J2( 2 d + k )- 2n - l ) p jQj+k + (2* - 2n - l)Q k , 

3 = 1 

that is, the equations \8J$ can be expressed as the following time- dependent Hamiltonian equations 

d Sk A^ = {h[ n) + 6h[ n) ,A M } + d c s a k n A (n *> 
with the Poisson bracket {, } defined by 

df dg dg df 



(86) {/,<?} = £ 



dP l dQ l dP, 8Q l 



and the derivative <9™ n is the explicit time derivative of Sk when the canonical coordinates Pi, Qi are fixed. 

Now let us deal with the tk -flows. Thanks to the fact that the coordinates u±, . . . , um exist and are such 
that l|77[) is satisfied, equations (|T5|) and l|76ll we see that the matrices MW correspond to the times (2i+l)tj. 
That is, the equations 

(2k + l)d? k A (n) = d x M {k) 

(2k + l)d tk A {n) = [M {k \A {n) ]+d x M {k \ k=l,...,n 

are compatible and hence we can express the higher order isomonodromic flows (|14(l as time-dependent 
Hamiltonian flows. 

We can now compute the Hamiltonians hP^ for the times t k by using Theorem 18 . 41 and (|76|l . 
Corollary 8.5. Let IC S and 7i t ~ be the following polynomials in A and j- respectively: 

n n— 1 

= E(4 n) + ^i n) )A 2fe , H t = + ^)(2A)" 2fe - 1 - i(H<-) + WW), 

fe=i fc=i 

where h^f 1 , 8hu^ are given by \85}) and respectively and the coefficients + 8H k are given by: 

(87) d\Tit = -2(2A)- 2 "((2A)- 2 /C s (l + T)- 1 ) + , 

with T is defined by \4b\j . Then the equations 

(2k + l)d tk A {n) - d x M {k] = [M {k \A {n] ] 

can be expressed as time- dependent Hamiltonian equations with the Poisson bracket given by f#6)) as follows: 

d tk A {n] = +5H^\A {n) ^ +d L t c k m A {n) , 

where the derivative <9™" is the explicit time t k derivative when the canonical coordinates Pi, Qi are fixed. 
In particular, the equation \1$ is Hamiltonian and it is equivalent to 

dQj _ dnj n) dPj _ dU k n ) 
dtk dP% ' dt k dQi 

where 



k = 1, . . . , n — 1, i = 1, 



H[ n) =Hl n) +6Hl n \ k = l,...,n-l. 
Proof. From J7HJ, we see that the Hamiltonians are related by 

n 

h ( k n) + 6h k n) = 4 k £ (2(* -k) + l)U(Ht\ + 6H™) A k t k (H^ + 5H™). 

i=k+l 
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This implies the following relation 

2K S = -(2A) 2 ((2A) 2n 9 A W t (l+T)), 
((2A) 2 "(2(2A)- 2 "- 2 /C s )) + = - {{2\) 2n d x H t {\ + T)) . 

Then, by using (HU, we obtain (JEZJl. □ 

Example 8.6. Let us demonstrate the results of this Section in the first non-trivial example: n = 2. We 
want to represent the t\ flow 

(88) id tl A {2) - 9 A M (1) = [M^\A {2) ] 

where 

4A 3 - 2w 2 \ -4wA 2 + 2w z \ - w zz + 2w 3N 

-4wA 2 - 2w z X - w zz + 2w 3 ~4X 3 + 2w 2 X 



16A 4 + 4A 2 (ti - 2w 2 ) - z - 2t lW 2 + 6w i + 2w 2 z - 4ww zz 
-l6X 3 w - 8X 2 w z + 4A(2w 3 - wt x - w zz ) - 2t x w z + \2w 2 w z ~ Sa. 

-16A 3 w + 8A 2 w z + 4A(2w 3 - wt x - w zz ) + 2t x w z - \2w 2 w z - 
-16A 4 - 4A 2 (<i - 2w 2 ) + z + 2t 1 w 2 - 6w 4 - 2w 2 + Aww zz 



as a time dependent Hamiltonian equation. 

We are now going to show that there exist coordinates iti, ... ,1x4 such that 

This gives the following equations 



(89a) 


3^4 2) 


= 39£ (-8w 2 + 4ti) = 12, 


(89b) 


3^>i 2) 


= 39£ (2(w 2 + 3w 4 - 2ww zz ) - 2t x w 2 - z) 


(89c) 


39 t >f 


= -483J> = 0, 


(89d) 




= 249 t >,=0, 


(89e) 


3^4 2) 


= 39£ {-Aw zz + 8w 3 - U lW ) = -8w, 


(89f) 


3dlbf ] 


= 39^ (2(w zzz - 6w 2 w z ) + 2t lWz ) = 2w z . 



We also need d^t x = 1 and d^z = 0. Of course, after imposing these two constraints, we have more equations 
then independent variables. (6 equations and 4 independent variables w, w z , w zz and w zzz ) We need to 
check that equations (|89a[l - lj89f|l are consistent. To see this, first note that equations (|89c|) and l|89d(l imply 



9 t > = 9 4 >, = 0. 
By substituting this into (|89a|l . we see that 

39£<4 2) = 39£ (-8w 2 + 4ti) = 12, 
which is tautologically true, so (|89a|) is consistent with (|89d|) and (|89c|) . By substituting i|89d(l and l|89c|l 



into (|89b(l . we see that 

39 t >^ 2) = 33£ (2(w 2 z + 3w 4 - 2ww zz ) - 2t x w 2 - z) = -2w 2 
dt 1 w zz = -\w. 



3 



Then by substituting these into (|89e|l . we find that 

3d^bf ] = 39£ (-Aw zz + 8w 3 - U x w) 
= -12d^w zz - 12w = ~8w, 

therefore l|89e|l is also consistent with the other equations. Now the last equation i|89f|) does not have any 
consistency issue and it would give us the following 

3&t 1 w zzz = -2w z . 

Therefore these equations are all consistent and u x , ... ,114 exist. 
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Using the results of Section El we can express the matrix in terms of the canonical coordinates 

-Pi) Pii Qi, Qi and using formula (|S5|) we get 

h (2) = -32zP 2 + 256P 1 2 P 2 - 256Pf - 2P 2 Q 1 Q 2 - P 1 P 2 Qj - 2a 2 {P 1 Q 2 - Q x ) - t 1 P 2 (128P 1 - 16tx) 

hf> = -32zP 1 + 256P? - 512P 1 P 2 - P 2 Q 2 2 + Q? + 2a 2 Q 2 - h (128P 2 - 128P 2 - 16Piii) , 
and using (|83|) we get for the shifts: 

Sh {2} = -P X Q 2 - Q u Sh (2) = -Q 2 , 

so that, using l|87|l . the Hamiltonian of the equation Pn'- 2 '' is 



n^ = h-h^ + Q 2 ). 

Id 

and the Hamiltonian of (|14Jl with k = 1 is 

M 2) -iQ(^ 2) + ^ 2) )-^i(4 2) + ^ 2 2) ) 

(2) 

For the sake of completeness, we can compute the shift in TL\ directly from the explicit dependence of 



Pi, P2, Qi, Q2 on 8%. From (|89a|l - H89f|l . we have 

d?Pi = ^(4 2) +4 2) )-^ 

ll 16 tlV 1 1 ' 24 12 48' 

d£Q2 = 0, 

dlQ x = -d^b™ -Q 2 dlP 1 = -9l. 

( 2 ) 

Therefore the we can find SH^ from the following partial differential equations 



This gives the shift 



sh[ 2) 



d5H {2) 


1 




dQi 


~4' 




d5H {2) 


Pi 




dQ 2 


12 " 


48' 


dSH {2) 


O2 




dPi 


12' 




d5H {2) 
dP 2 


= 0. 




i_ Qi 


P1Q2 


hQ 2 


4 


12 


48 



which in fact agrees with what we obtained above. 

Appendix A. From the mKdV Lax pair to the isomonodromic problem 
Let us consider the Lax pair of the mKdV hierarchy: 

— = 7V$, = 7V fc $, 

ox ot k+1 

where 

N _ { -( v \ n - ( T.f=i A f ] v v;>/- \ 
{ v C J' Nk ~{ e 2 L 4 ] t j - eS 1 4 } v ) ' 

where the coefficients A^ k \ B^ k \ are given by the formulae i|16|) with Ci replaced by IZi, z replaced by 

'a 



(k) 

x and w replaced by v. Observe that in N k the terms B a are nonzero. They are 



B { k) = -{d x + 2v)TZ k . 



Details can be found in |10j . 
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2!) 



By the self-similarity reduction (jl 1|> , and defining 

A 



we obtain: 



(2n + l)T nH 



[(2n + l)T n+ i 



dz 



{2k + l)t k - + X-, 



i 



9T fe 



- t k 



Now, defining 

we get 
(90a) 

(90b) 
(90c) 



B=[{2n + l)T n+1 ] — N, 
dz 

a* 

~d\ 



dx dz 
M« = [(2n + f)T n+1 ]^TT7V fc , 



fe= 1,. 



6*, 



To show that l|90[l coincide with (|f 511 first observe that 





(91) 



M k = M (fe) - ( 



{d z + 2w)C k 



(8 Z + 2w)C k 




as it follows from the definition of M k and Nk- Moreover 
(92) A {n) = M {n) + s ^t l M {l) 



zB, 



1=1 



because 



-{d z + 2w)C n 



n-1 

£ 

i 



ti(d z + 2w)Ci + zw = —a n 



thanks to JT2J). 

It is now clear how to prove that the equations (|1 T|) . I|18|l and (|I9|I are indeed consistent. This follows 
from the fact that the mKdV flows commute, so that 

ON dN k 



and 



dT k+1 
dN t 



dx 
dN k 



N k = [Nk, N], 
= [N k ,Ni\. 



8T k+1 dT l+1 

The equations (|17|l . (|18|l and l|19|l then follow automatically from the relations i|92|) . I|91|l . 
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